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To avoid problems with infinite measure, the functional integral for harmonic oscillator can be cal- 
culated by time - slicing method with continuum limit procedure proposed Gelfand and Yaglom. In 
previous article we proved by nonperturbative calculation the generalized Gelfand- Yaglom equation 
for anharmonic oscillator with positive or negative mass term. In this article we prove by step- 
by-step the calculation of the correction function to the Gelfand- Yaglom equation for an-harmonic 
oscillator. 



Introduction 



Let us demonstrate the meaning of continuum limit procedure proposed by Gel'fand- Yaglom on the example of the 
harmonic oscillator. In Euclidean variant of the theory the continuum functional integral for harmonic oscillator can 
be read: 



Z = J [Vip(x)} exp(-S) 



Where the euclidean action is: 



S = J d r 
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In time-slicing approximation of Wiener unconditional measure functional integral we calculate the N— dimensional 
integral [l[: 
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where A = P/N, and b, c are the parameters of the model. The unconditional measure N dimensional approximation 
(fixed ipo and integration over ifN ) can be evaluated explicitly: 
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where LUi is defined by recursion 
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with the first term 

uj = 1/2 

where 

A = 



bA 2 /c 



2(l + 6A 2 /c) 



1 



2(l + 6A 2 /c) 



The recurrence relation for functional integral is introduced by time - slicing procedure. Explicitly this recurrence 
is represented by factor un. Following this recurrence, we can prove the difference equation for inverse square root of 
the finite dimensional integral. The continuum limit of this difference equation give the non-trivial value for functional 
integral. Based on this procedure the continuum limit N — > oo of the Z® N can be defined by: 

Z°((3) 



It was shown by Gelfand and Yaglom that F(j3) is the solution of the equation: 

8 2 2b 

jr( T ) = _ F( T ) , re (0,0). (4) 

For the harmonic oscillator the equation (j4|) can be calculated by taking the continuum limit of the difference 
equation extracted from the recurrence relation for Z^- 
The general solution of the above equation is: 



F(t) = d cosh [J— t + C 2 sinh [J— r (5) 





where C\,C2 are the constants fixed by initial conditions of the unconditional (when ipo is fixed and ip^ is free), 
or conditional (when both (fo,<fiN ar e fixed) measure functional integral. From the analytical form of the result for 
functional integral follows the clear interpretations of the calculated quantities as energy levels and others. 

In this article we report on the attempt to evaluate non-perturbatively the functional integral for an-harmonic 
oscillator with positive as well as negative mass squared term. Our aim is to calculate the N— dimensional integral 
by another method as well-known conventional perturbative calculation. We find the result suitable to evaluate by 
the recurrence procedure the difference equation. Following the idea of Gelfand - Yaglom we define the differential, 
Gelfand - Yaglom type equation for the quantity y(r) related (similarly as F(r)) to the unconditional measure 
functional integral. The differential equations for y(r) reads: 



8 2 8 8 (lb 8 2 8 \ 

^(r)+4- y (r) _l n 5(r)= 2/ (r) {- - 2— 2 In S(r) 4(- In S(r)) 2 J , (6) 
We shall evaluate the function S(t) analytically in this article. 

As we explain in the text, our result possesses the form of an asymptotic expansion power series. Nevertheless we 
evaluate precisely all finite difference mathematical objects. In the article quoted as "article I" and in this article 
we explain and prove all analytical calculations. 

For an-harmonic oscillator the correction term is —2^- In S(t) — 4( In S(t)) 2 to the equation for the harmonic 
oscillator. It would be desirable to recover it using more general arguments and to see its form for some more general 
classes of potentials. 

This article is organized as follows. In the next section we resume the article I. In the third section we evaluate 
the function S (r) by recurrence relation from the result of N— dimensional integral. In the fourth section we discuss 
some preliminary conclusions of our calculation. 
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Resume of the article I. 



Our aim is to solve the problem of evaluation of the continuum unconditional measure Wiener functional integral: 

Z = J [D<p(x)]exp(S) , 
where continuum action possesses the fourth order term: 



s = J iT 





c/2 



( Mr) 

V dr 



bip(T) 2 + aip(r) 4 



The functional integral Z is defined by limiting procedure of the finite dimensional integral Z^: 



(7) 
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where A = /?/iV. Then, the continuum Wiener unconditional measure functional integral is defined by the formal 
limit: 

Z = lim Z^v • 

The first important task is to calculate the one dimensional integral 



Ii = dx cxp{ — (ax A + (3x 2 + jx)} 



(9) 



for Re a > 0. 

Standard perturbative procedure rely on Taylor's decomposition of exp(— ax 4 ) term with consecutive replacements 
of the integration and summation order. The integrals can be calculated, but the sum is divergent. 
We propose the power expansion in 7: 



+00 



Ji= y*(_2T / dx x n exp{-(ax A + f3x 2 )} . 



n=0 



The integral is given in terms of the parabolic cylinder functions: 



(10) 



D- m -l/2(z) = 

The integral I\ then can be read: 



-z 2 /4 poo 



where 



and we have used the abbreviation: 



dx x m 1/2 exp{--x 2 - zx} 
r(m+l/2)7 yx 2 S 



h = 



r(i/2) ^ c 



v ' m=0 
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It was shown, that sum in Eq. (fT2|) is convergent and for finite values of the parameters of the model this sum converges 
uniformly. 

Applying this idea of integration on the N dimensional integral (JSj) integral we have the result: 



Z N — 



N 



Y[2(l + bA 2 /c)uj % 



i=0 
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with 



Sn = 



e n 



(p) 



2k, 



{2ki) 



■Tih-i+ki + 1/2) y/UTiV. 



-fe«-i 



/2 (*) 



(13) 



(14) 



where the constants and symbols in the above relation are connected to the constants of the model by the relations: 
fco = k N = 0, p = (1 + bA 2 /c)-\ z = c(l + 6A 2 /c)/v/2aA3, w, = 1 - A a /wi-i, ^o = 1/2 + 6A 2 /c, A - 2{1+b \y c y 
In the formula for Sn is useful to mention that: 

- p is independent of the coupling constant; 

- only the argument z of parabolic cylinder function is coupling constant dependent; 

It was shown, that for finite values of the parameters of the model and one summation index fcj — > oo the fcj — i/i 
term of the sum approaches to zero as 



fei! exp(vfci) 

where a and (3 are finite numbers. This asymptotic is sufficient for a proof of the uniform convergence of the series 
for Sn not only for single fej, but for arbitrary tuple {ki} of indices. 

Following the idea of Gelfand and Yaglom the functional integral in the continuum limit is defined by the formal 
limit 

lim Zjj- = — ■ , 



where F((3), for our case, is the solution of the differential equation: 



dr 2 



F(t) 



4 1 i F{T) 



d_ 



]mS(t))=F(t) 



2b 



2^1n5(r)-4(A lniS(T))2 



re (0,0) (15) 



calculated at the point j3 (the upper limit of the time interval in the continuum action) with the initial conditions: 
F(0) = 1 and <9F(t)/i9t| t=0 = 0. Evaluating the continuum limit of the difference equation we use the convention 
for definition of the continuum variable r = n A . 

The function S(t) is given as the continuum limit of the Eq. ([M]) 

S(t) = lim Sn ■ 



S(t) 2 



Equation (| 1 5[) can be simplified by the substitution: 

F(r) 

For y(r) we find the equation: 
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_y( T ) = y( T ) ( _ 



2b 



(16) 



accompanied by the initial conditions: 



y(0) = 5(0) 2 



dy(r) 



Ot 



T = 



5 



For harmonic oscillator we have 5 = 1. 

To evaluate Sn, we must solve the problem how to sum up the product of two parabolic cylinder functions in Eq. 
(|14p . The parabolic cylinder functions are the representation of the group of the upper triangular matrices, so we 
implicitly expect the simplification of the product due to a group principles. This problem was not solved completely 
yet. We adopt less complex method of summation, namely we use the asymptotic expansion one of parabolic cylinder 
function, with precise sum over ki of the rest of the relation, containing the other function. In the algebra such 
relation for precise summation is available. Surely, the result is degraded to the form of an asymptotic expansion only, 
but still we shall have an analytical solution of the problem. This procedure is widely discussed in article I, here we 
repeat the result: 

Z N = m [2(l + 6A 2 /cK]| S N , (17) 
*A =E7T7S^A3,(A)^, (18) 



where the symbols (A) 2 "* satisfy the following recurrence relation: 



= £ ( X ^ £ (db) < A - °-" + " 

A=0 V 7 W A-1 i=max [0, 2A-p] A z A 1 / 

The recurrence procedure begins from: 

We repeat also the definition of the value z, where the dependence on the coupling constant is hidden: 

_ c(l + bA 2 /c) 

As follows in the calculation the important role plays the objects tVi defined by recurrence as 

A 2 

Wi = 1 , 

OJi-l 

with the first term for unconditional measure integral: 

6A 2 /c 



(19) 



w = 1/2 + 



2(l + 6A 2 /c)' 

For the forthcoming calculation we will use the more convenient variables, introduced in Appendix 2 of article I: 

Qi = w 1 x l + w 2 y\ 



Qi = w\x l ~ wiy\ 



where 



^2 i 4 + V3i" 1 '^^"Vii- 1 ' Xy=1 ' 



2B 2B 

w i = 1 H ; ! w 2 = 1 



VI - 4A 2 VI - 4A 2 
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1 bA 2 /c 
A= n/ , , tA ,„ , B = 



2(l + 6A 2 /c)' 2(l + 6A 2 /c)' 

We also use the fruitful identity following from the definition of the n — th convergent of the continued fraction, 
calculated in article I: 



A Q 



fc-i 



Wk-i Qk 

The symbols a\ were defined by need to rewrite the Pochhammer symbols (k + 1/2) . in the form: 

min (j,k) 

(fc + 1/2) = a{k(k-l)...(k-i + l) , 



i=0 



and by help the recurrence procedure we found: 



* J (1/2),- ' 



In this article we explicitly evaluate the recurrence relation for (A) 2 ^ and we calculate the continuum limit of the 
function 5a- 



Evaluation of the recurrence relation 



We rewrite the recurrence relation p9[) into more convenient form for consecutive calculation. We introduce the 
quantities Qk by the identity: 



.4 



!k-l 



Uk-i Qk 

We replace the summation index i by the summation index j defined by: 

i = 2A - j . 

Finally, we interchange the order of summations over indexes j and A. We read: 
{AQ K Q K - X f d -P{k)l d d _ v = 

- E (AQ a qZi)^ E (^a- 2 Qa-i) 2A ^ ( A - ( A ) ^l-i)^ • P e < 0. M > ( 2 °) 

The right hand side of the equation is (2d, p)-th matrix element of the products of three matrices. For fixed d and 
p, on the left hand side of the equation, we read only the d - th column of a matrix, which is recurrently tied to the 
matrix in the center of the product on left hand side. We can use the notation: 

A ) = E E ^(A-l)C^ iA (A-l)M^ iM (A-l). 

J=0 A=[2±l] 

The definition of the matrices A d (A - 1), C d (A - 1), M rf (A - 1) is the following: 

1. The A d is the lower triangular matrix with the zeros over the main diagonal of the dimension (2/j, + l)(2/i + 1). 
The principal minor of the dimension (2c? + l)(2d + 1) is non-zero only with the elements: 

a 2d-p 



|A d (A-l)| . =— 

1 v >j m {AQ\Qa-i) p 
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2. The C d X (A) is the upper triangular matrix with the zeros under the main diagonal of the dimension (2/z+l)(/i+l). 
The nonzero elements form the main minor of the dimension (2d + l)(d + 1) with X th column: 

C& A (A) = (AQaQa-^-v (A)l x x _ p , 

where p = 0, 1, 2A and A = 0, 1, d . 

3. The Mf k is the upper triangular matrix with the zeros under the main diagonal of the dimension 1)(/Lt+ 1). 
The nonzero elements form the main minor of the dimension (d + l)(d + 1) : 



;, fc (A-l)=(*)(Qi_ 1 ) fc - A , d>fc>A>0. 



To evaluate the matrix C(A), we must calculate X(A) for all dimensions up to fi, for each dimension to extract the 
last column of matrix X(A) and from these columns to compose the matrix C(A). We define such linear operation as 
follows: 

1. Let A d and M d are the matrices of the dimensions (2/i + l)(2/x + 1) and {fi + 1) * fi + 1) respectively. Let C d 
is the matrix of dimensions (2/i + + 1). This matrices possesses the nonzero main minors of the dimensions 
(2d+ l)(2d+ 1), {d+l)(d+ 1), and (2d + + 1) respectively. 

2. Let supM is the supermatrix possessing on the d — th place of the main diagonal the matrix M d the same is 
defined for supA . 

3. The supC is the supermatrix with d — th diagonal element of the form: 

d 

C d (A) = E X d ~* A (A) . 

*A=0 

4. Matrix X d is the one column matrix defined by the relation: 

X d (A) = X d (A) * P d , 

where P d is the projector of the d-th column of the matrix X d (A) into d-th column of the matrix X d (A). F d is matrix 
with only nonzero term 

5. The matrix X d (A) is defined by relation: 

X d (A) = A d (A - f ) * C d (A - 1) * M d (A - 1). 

6. Then, for C d (A) we have the result: 

d 

C d (A) = E A d_iA (A- 1) *C d - lA (A- 1) *M d ~ 4A (A- 1). 

i A =0 

7. After evaluation of the full recurrence we find: 

d d—i\ d — iA — <a-i i;t 

c d (A)-E E ••• E ( 21 ) 

i A =0iA-i=0 i2=0 

{A d " lA (A - 1) * A^^-^-^A - 2) * • • • * A d -»A-i A -l i2(]_)J * C d-»A-»A-1 »2( 1 ) * 

|j^d-»A-»A-i » 2 ( 1 ) * . . . ♦M^^-^-^A- 2) * M d_ * A (A - 1)| 

To evaluate the product of two consecutive matrices from multi-product: 

{A d ~ 4A (A - 1) * A^-^-^-^A - 2) * • • • * ^d-iA-iA-i ^2^)1 

we use the two identities for the summation over the index j: 



2/3-, . 2/ 2 -a _ 2 - 2( p-A) (4J 2 -2A)! / p - A \ (4J 3 -2j)! 
«2/ 3 -P «2/ a -j - (4/3 _ 2p)!(p - A)! V j - A J (47 2 - 2j)\ 



and 



(4/ 3 -2j)! 



a 4J 3 -4/ 2 ( e «s-2j)| e=1 



(4/ 2 -2j)! 
We introduced the abbreviation 

Ij =d - (ij + H h u) 

Then, for product of two lover-triangular matrices we have: 



2/ 2 -A 
2/2 -3 



)i-A 



2 -2(p-A) 



(4/; 



(4^2 -2A)! ^.L^f/p-AW e 2 
,-2p)!(p-A)! e | l j J tlU'- A ; UO3Q2 



Q2J \AQ 2 Q 1 



In the above relation the summation over index j can be performed explicitly. By the substitution: 
we find: 



^{A^( 2 )} j)i .{A^(l)}. A = 2- 



2(p-A)_ 



(4/ 2 -2 A)! 



-2 4i2 D i2 



1 



1 



(4/ 3 - 2p)!(p - A)! « I C p - 2/3 \ AQ3Q2 AQ 2 Q! 
where is the differential operator calculated from <9^ given as: 



p-A 



Df = 3/4<9f + 3^6»| + £ 2 <9| 
For the resulting product of all matrices A / (fe) we find: 

{A d - lA (A - 1) * A^-^-^A - 2) * • • • * A d -»A-»A-i ^ = 



2 -2(p-A) (4/ 2 -2A)! q4( I2 +» 3 + ---+»a) - 

(4J A -2p)!(p- A)! 



< m=2 



^p- 2 / ro+1 ^q 2 Qi 1 aq 3 q 2 

Evaluating the product of matrices: 



1 + 6 6' • 'Ca-i 



p-A 



(a«i£ m ->l) 



»A— 8A-1 12 



(I)*---*! 



- 1 (A-2)*M d - 4A (A-l)| 



we use that: 

M /j+1 (j) is one-column matrix with I 3 +\ non-zero elements in j + 1 column: 
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_ [ Ij+i j g 4 ( / j+ 1 " A ) 



where A = 0, 1, • • • , Ij+x- Product of such matrices is one-column matrix with the elements: 



-IA— «A-1 «2 



(1) 



nd— ia— tA- 



1 (A-2)*M d -' 1A (A-l)| 



A, /a 



/ 2 W /a 
A J U 



\ n 4(/ 2 -A) n 4(/a-/2) n 4(I A -/ A _i) 
J | Vl ^2 '"Wa-1 



From definition (|19p of the recurrence steps we have for the matrix C /2 (l) the nonzero elements: 



{C*(l)} 



n 4A 



with the conditions for indices: 



< j < 2A < 2I 2 < 2fi 



Collecting all partial results together, inserting them into Eq. ([2Tj) and remember that for function 5a defined in 
Eq. (fT5|) only matrix elements {C(A) 2m } 2 ^ 2 ^ are important, we find the result: 



{C(A) 2 ^ 1 } 



2fi,2ft 53 



H /i— !*a M— «A~ 'A-l- 

E ••• E 

lA=0 i A -l=0 42=0 



7 A \ n 4(/ 3 -/2) n 4(/ A -J A _i) 



^ * (4J 2 - 2j )!2 4 fa+-^) [A 

> 2-, (On _ 1\\ 0(4u-2j) 1 11 ^€ 



J=0 A= |2+i| V P ^ 



v m=2 



1 



2^-2I„ 



1 



£2 • ■ • £ 



A-l 



AQaQ 



A-l 



(aHU-1) 



1 



^ 2 1 „ 2A n 4A n 4 ( /2_A ) 

\ I a 2A-j Vl 



AQ Q 1 

The asymptotic function 5a can be read: 



(23) 



J 



5 A = £ 



fi\ (2z 2 A 3 y 



2fj,,2p. 



(24) 



because {C(A)^} 2 ^ = (A) 2 ^ 

Due to the analytic form for {C(A) 2m } 2 ^ 2 ^ we can express 5a in the continuum A — * limit as well as in asymptotic 
(i — > 00 limit. 



The continuum limit 



In continuum limit we must take into account that: 

- in Eq. the fi - th term {C(A) 2 ^} 2 ^ 2 ^ is multiplied by A 3 ^. 

- the sum 



1 



AQ 2 Q 1 AQ 3 Q 2 



£2 • ■ ■ £a-i 
AQaQa-iJ {aUU ^ 1} 
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is given by exact formula: 



1 (Qa_Qo 
2Aw!W 2 (x - y) \Q A Qo 



From the relation 



2Awiw 2 (x - y) y/2b/c A' 



we deduce that to the leading term of Eq. (|23[) contribute the terms with summation index j — with the contribution 
proportional to I I 

- to obtain the additional necessary factor (^^j the leading term must be composed from the contributions where 

fj, summation indices of the {i 2l ia, • ■ ■ , *A-i} are equal to 1. Let it be the combination, say, {i ni , in 2 j ■ ' " ? *n }• The 
sum of these indices is equal to fi, and in such case also I 2 — 0, because of definition I 2 = /i — i ni — i n2 — ■ ■ ■ — 
The contribution to the body of principal formula (I23|) then can be read: 



'»2 

n i 



I rip 



) 4 



(25) 



(2/x)! 



II 2 M -2/„ m + 1 



s m.—l 



l 



6 



aq 2 Qi aq 3 q 2 



6 • • • £ 



A-l 



A-l 



2/' 



(aZ2£-»l) 



We must sum over all such combinations, this can be done by /i summations. Every summation is proportional to 
A -1 . Therefore leading term proportional to A~ 3m can be achieved only if /i different indices i n = 1. Taking into 
account, that the difference of two consecutive I ni is one, we can rewrite the dominant contribution in the continuum 
limit into the form: 



{C(A) 2 "} 



A-£i A-/J+1 A 

E E ••• E 

«i=2 ri 2 =rii + l 78^=71^-1+1 



(26) 



s m— 1 



«»m 2 M -2/„ m+1 
fern 



1 



6 



AQ 2 Qi 



6 ' ' ' £a-1 



A-l 



2/' 



(aZZ^-»l) 



The effect of the operator 

Let us evaluate the operation of the operator D x accompanied by operator's variable term: 



1 



lim D x —,. 

Xi -+1 X i 2(z 3 + i+---+ lA ) 



(27) 



ij times on the n — th power of the function 



fj = dj-i + {x 2 - ■ -Xj)bj 

where: 
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Qj+iQj QaQa - 1 

This definition reflect the fact, that before application of D Xj we provide the limits Xj+i — * l,---,X\ — » 1. The 
following limits will be used: 

lim (x 2 ■ ■ ■ Xj) = (x 2 ■ ■ ■ Xj-i) 

Xj — >1 

lim fj = = aj_ 2 + (x 2 ■ ■ ■ Xj-i)bj-i 

Xj — >1 

We meet very important feature of the application of the operator D x . 

1. Applying D x the first time for first nonzero ij, all ij+i, ■ ■ ■ ,i\ are zero, then factor 2(i +... +iA) = 1. 

Xj 

2. We find in the evaluation, that there always appears the terms, killing the variables in denominator of (|27|) in the 
next steps of calculation. 

From the practical reasons, our aim is to express the resulting formula in the form where the dependence on the 
next derivative variable, Xj-± is in the function fj-i only. Wc find: 

}^ D % ^J + ... +iA) fr= (^) 

MIN , 

[\x 2 ---x j - 1 )b j \" : > 2^ I 

1=0 ^ 

where 



[(x2-~x j -i)b J ] 2i * J2 I a 0" * '.V.I.Ul.X-.n.iy.b, h, L ) 



J(i ,M» ; M,v,-^i:(r/)( 2 ^ 2 )(S^(a" 

p — I 

MIN = min (2i J ,n - 2y). 

The proof of this formula is given in Appendix A. 

We see, that in the term \{x 2 ■ ■ ■ Xj_i)bj] 2l i is the variable power xj" 3 canceling the same power of the variable 
in the denominator of the operator effecting over variable Xk, to left the evaluation simpler. It can be shown, that 
function J is proportional to the Gegenbauer orthogonal polynomial following the relation for the Jacobi orthogonal 
polynomial P^'^ (see e.g. Prudnikov Q): 



k 

k ) \n-k ) \x-l) ~~ (x - 1}' 

If the /3 is a half number, the Jacobi polynomial can be expressed by help of the Gegenbauer polynomial. 
The result of two consecutive application of the operator D x . and D x _ t is: 



MINI , s MIN2 , s 

[{x2---x j ^)b j ^r^[{x 2 --- Xj ^)b j ]^ 7 m' E V) (o ''- 3) ' 

(f^yi-Ki-Ki-i-t-i J(Z,M/JVl;n,ij;6j/6j_i) J(l,MIN2;n-2i j - Mj-U&i-i/fy-a) 



MINI = min (2ij,n - 2ij), 



MIN2 = min (2^_i, n - 2ij - 2ij_x - I). 
On the result of application of three operations the nonlinear character of the our result is clearly visible: 

\ lim D^rl ,,. \ -. , 1 I Um D%r\ ,,. j , 1 J lim Z3*?'. „.. * - , 1 (/,•)" = 
[(a* • • • Xj-dbi-*]*'-' [(x 2 ■ ■ ■ x^bj.^-i [(x 2 ■ ■ ■ x^bj] 21 * 

MINI /„. x /■ \m MINI , v , , s i-m MINI , . \ 



, . ,(';.-.■■> / ' — ' \ l — m i \ n. ■> 

m=0 x / \ J z / i- m \ \ ./ 



MIN3 , 

2l 



2=0 



where 



E T 2 M"(/,-3r 2 ^ 2l - l2! - 2 " ! ^ JO.MJJVl;,,^;^/^!) 
J (7, MIN2; n — 2ij — I, ij-x; bj-i/bj-2)3(v, MIN3; n — 2ij — 2ij^\ — m — I, ij-2', bj-z/bjs) 

MIN3 = min (2ij_2, n — 2ij — 2?j_i — 2ij~2 — m — I). 

The above relations give first three terms of the asymptotic expansion for function S\ ([24|) . We find: 
For fi — 1 only one summation index ij is nonzero, then 

I A 

{C(A) 2 } 2 2 = J±- E Qt bt 3(0, 0; 2, 1; 6 fc /6 fc _ x ) 

^ k=2 

For [i = 2 two summation index ij are nonzero, then 

. A-l A 

{C(A) 4 } 4 4 = T^Tj E E Q^p^p J(0,2;4 ) l;6 fc /6 p )J(0,0;2,l;b p /6 p _ 1 ) 

^ fc=2 p=k+l 

For = 3 three summation index ij are nonzero, then 

^ fc=2 p=fc+l q=p+l 

E ( / ) ( X - ^) J ( Z ' 2 5 6 - ^ bk ^ J (°' 2 - i; 4 - 1, 1; J(0, 0; 2, 1; VVO 

The continuum limit is introduced by prescription: 

Afc — * x, AA^t, 
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Fig. 1: b dependence of the continuum function S(a, b, c, r) for fixed values a, c, r. The first three nontrivial terms of the 
asymptotic series (|24[) were used. 



A 1 ,T 



A 

k=2 

When A — » AT, then r — » /3, where /3 is the the constant of the model, and N = 
In continuum limit we obtain: 

Qfe — > 2 cosh(7x) 



6fe — > — — (tanh(7r) — tanh( 7 x)) 

A7 



where 7 = y / 2b/c, b and c are the parameters of the model. The continuum limit of the relation (|24|) we will call 
S(a, 6, c, r). To ilustrate the analytical form of the result, we show the first nontrivial term (ju = 1) : 



{C 2 (a,b,c,T)} 2 2 = ^3 [3 7 Ttanh 2 ( 7 T) + tanh( 7 r) - 7 t] (33) 

For the higher we have the analytical formulas also as the results of algebraic evaluation by Mathematica. The 
continuum function S(a, b, c, r) for the first three nontrivial contributions is shown in the Fig. 1. 

The corresponding term for the Gelfand-Yaglom equation, — 2<9 2 ln(5(a, b, c, r)) — 4(<9 M ln(5(a, b, c, t))) 2 is shown in 
the Fig. 2. 



The leading divergent term in the limit n — > 00 

In this limit the terms ([23)1 are divergent. We are going to evaluate the leading divergent term for each /i. To 
provide this, we must to evaluate the double sums in Eq. 



2/2 



E E 



(4/ 2 -2j)! 2 4 ( l2 +- +lA ) 



^ {2a-i)\ 2^-23) 



1 



2ft— 2J„ 



^Q 2 Qi 



£2 • ■ • £a-j 



A-l 



(? m -l) 



^4<3o<3i 



^ 2 1 ^ 2A Q 4A Q 4 ^ /2_A ' 1 



(34) 



The details of the calculation are explained in the Appendix B, we obtained the result: 
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tau 



Fig. 2: r dependence of the continuum function — 2c? 2 , ln(S(a, b, c, r)) — 4(9 M \n(S{a, b, c, r))) 2 for fixed a,b,c. The first three 
nontrivial terms of the asymptotic series (|24|l were used. 



{W 2 %, 2M = (1/2) 2M (- 



1 



1 



2/i 



«j 3 + q ;)^2mi/2) 2 „(^) 



2,, 



!2 / ,2„ '^V^Qi AQ A Q A I 

For the leading divergent term — > cx) of the asymptotic series for 5a (I24[) we finally have: 



2/i 



(35) 



The series of this form is an asymptotic expansion of the parabolic cylinder function of the index —1/2 and the 
argument 



-l A ( tanh( 7 r) 



V c 7 



Gelfand - Yaglom equation for anharmonic oscillator with mass zero 

As a test of our calculation we evaluate the energy levels of anharmonic oscillator with zero mass. The Gelfand 
Yaglom equation for this case will be obtained by limit b — > 0, and we have: 



y"(T)+4go(T)y'(T) = fo(T)y(r) 
For the functions /o(t) and go (r) we evaluate from Eq. (|24p for first three non-trivial terms: 



(36) 



/o(r) 



9o(t) 



dr 2 



In 1 



ax 3 43 /ar 3 \ 2 8111 /ar 3N 3 



30 



1890 



d ( ar 3 43 /ar 3 \ 2 8111 /ar 3 ' 3 
rf7 n I 1 + 30 (~ ~ 1890 l^2 



d I ar 3 43 / ar 
— In 1 — 



3 .1Q /_3\ * glll / ar 3N 3 



V c 2 30 V c 2 / 1890 
For first order in the coupling constant a, we have the equation: 



(37) 
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„ 3.4ar 2 3.4ar 
y 0) -2— y (r) = — — y(r) 



(38) 



This equation can be solved analytically, (see Kamke @, Eq. 2.60). By substitution y(r) = u(t) exp (2ar 3 /c 2 ) we 
find the equation: 



u"(r) 



T 4 m(t) 



The general solution is expressed as the linear combination of the Bessel functions: 



,y T ,■ 2aT3 ^ /■ 2aT3 ^^ /2ar 3 \ 

»(r) = V? ( CiJr/e^-^-) + C 2 r 1/6 (z— — ) 1 exp ( — ^ J 



The constants Ci and C2 will be fixed from boundary conditions. For z~0we follow the identities 

(-1)' 



/ z\ u \ ( — l) m /z\ 2m 
= (2J ^ m! T{m + u+l) \2) 

m— 



and 



Y v (z) = — — - [J„(z) cos(^tt) - J-^)] 
sin (vn) 



(39) 



Inserting to the equations y(0) — 1 and y'(0) = we find: 



,\ 1/6 



C i = - -5 T(l - 1/6) sin (tt/6) 



1/6 

C 2 = I =1 ] r(l - 1/6) cos (tt/6) 



For evaluation of the energy of the ground state we need y(r) for r — > 00 . For this limit we can use the relations 



Tf\ 2 ( 2V + 1 
J V [Z) = \ COS Z 7T 

V TTZ \ 4 



Y v (z) = \ J — sin z — 



2v + l 



and finally we have for the functional integral in this limit: 



1/6 

3 ) r (! ^ V6) 



1/2 



sm 



/2ia/? 3 7T 



g I exp 



2a/3 s 



1-1/2 



The unconditional measure functional integral is the partition function for the model solved, for the harmonic 
oscillator there is the simple relation for the energy of the ground state: 

Eq = - lim (\\nF{!3) 

/3->oo \p 

By the direct application of this relation we obtain zero for the ground state energy of the an-harmonic oscillator with 
zero mass. We find the nonzero result for the definition: 
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Conclusions 



In this article we calculated step-by-step the correlation function for differential Gelfand-Yaglom equation an- 
harmonic oscillator with positive or negative mass squared term. We stress, that generalized Gelfand-Yaglom equation 
is the non-perturbative equation, the correlation function is evaluated in the form of the asymptotic series. The 
analytical form of the correlation function enables us to evaluate the continuum limit as well as asymptotic limit. The 
continuum limit can be used for evaluation such physical quantities as energy levels. The asymptotic limit cam be 
used for attempts to sum the series by Borel's method. 
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Appendix A 



We are going to apply the derivative operator 

D Xj =3/2d*. + 3x j d*.+x 2 j di j (A.40) 
i times on the n — th power of the function 

fj = Oj-i + O2 • --x^bj 

where: 

a - 1 I X2 \ x 2 ---Xj-i 

3 Q2Q1 Q3Q2 QjQj-i 
h 1 1 



Qj+iQj QaQa - 1 

We see, that bj is a constant from point of derivative operator D Xj and Oj_i also. Therefore for (D Xj Y /" we find: 

(D Xj y /; = 

min(2i.n — 2i) , , . 

{{x^.-x^r y, c-Wm) °r-i (n _ 2i- mv {n - 2i + 1/2) *-™ f ?~ 2i ~ m (A - 41) 

m=0 v / ^ ; ' 

The above result is nonzero only if 2i < n. The variable Xj—\ is in the term fj and dj-i as well. By the identity 

ctj-i = fj - (x 2 x 3 ■ ■ ■ Xj-i)bj) 

we introduce summation due to binomial expansion of Oj_i into Eq. (|A.41[) over index p = 0, 1, • • • , m. We exchange 
the order of summations and apply the identity: 

2i\ ( m \ _ ( 2i\ ( 21- p 
m J \ p J \ p J \ m — p 

Finally, performing the limit Xj — ► 1 we find: 

min(2i : n — 2i) . > 

(D X] Y f?\ Xj=1 = {{x 2 x^--x j ^)b j f i ( p ) {{x^---x j . x )b j )Yf^ i - p 

rni7i(2i,n — 2i) > > ^ 

Y (-if-™ ( 21 ~P ) £ ( n _2i+l/2) 9 . (A.42) 

m=p \ s \ 
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Now, we can perform the summation over index m. For the case min(2i, n — 2i) — 2i we use the identities: 



71 ' 

= lim dl* + ™ (x n ) (A.43) 

(n — 2i — my. 

and 

(n - 2i + 1/2), . = lim df~ m K > l0 (A.44) 

Inserting this into sum over m in Eq (|A.42[) . replacing the order of limit and summation, as well as the summation 
index m by I = m — p we have: 



When rmn(2i, n — 2i) = n — 2i we find the identical result. We proved that: 

(D Xj y /;u, =1 = 

min(2i,n—2i) / ■ \ I 

((x 2 z 3 • • • ^ ( * J ((x 2 x 3 • • • x^b^f^^ (n _^_ p) , (P + l/2) 2l _ p (A.46) 

In the sum over index p we recognize the formula for the Jacobi orthogonal polynomial. By help of the identity: 

/ 2% \ n\ = (2Qln! / n - 2i \ ( 2i - 1/2 \ 

\P J (n-2i- P y. [P 1 >2 *-p (n-2i)\\ P ) \ 2i - p ) 

for the case min(2i,n — 2i) = 2i we find the similarity with identity (see e.g.Prudnikov 0]): 

t("; n )(^)(^)'-^^' 

Where {x) is Jacobi orthogonal polynomial. 

The expression in Eq. (|A.46[) we simplify further for application of the operator (D Xj l ). In Eq. (|A.46[) the variable 
Xj-i is in two terms, therefore we simplify that equation by the identities: 

(x 2 x 3 ■ ■■x j ~i)b j = -p— (x 2 x 3 ■ --Xj^bj-i, 

and 

(x 2 x 3 ■ ■ ■ Xj-^bj-x = fj_i - a^ 2 . 
Now, the variable Xj-i will be in fj-i only and for application of the operator (Dx^^ we prepare the relation: 



min(^n-2i) p p n!(p+l/2)_. 



v* AW ^: V " / (n-2*-p)! 
By exchange of the order of summation we finally red: 

min{2i,n-2i) . . mm(2i,n-2i) . . . . p | / , -, /p\ 

( ( , 2 X3 • • • E ( ; ) E ^ " ^ ^ 2l ~ p 



V~ v ) \b 3 -ij (n-2i-p)\ 
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In this form is result of application of operator (D x .) 1 suitable for the next evaluations. 



Appendix B 



We are going to evaluate the leading divergent term. To provide this, we must to evaluate the double sum in Eq. 



HMD: 



^ ^ (4/ 2 -2j)! 2^+'3 + -+i A ) f a 
^ ^ (2u-j)\ 2( 4 ^-2j) I 11 «" 



1 , 



1 J{(t)^Q^- A) 



To accomplish this calculation we change the order of the summations and we apply the identities: 



(A.48) 



(4/ 2 - 2j)\ = 2< 4I »- 2 fl (27 a -i)! (l/2) a . 



2/2 -j 



,2A 



2A \ (1/2) 



2X - J (V2) 2A _ 



2A 



(1/2) 



2/j-j _ q2I 2 -2\ 



i 



2A-J+1/2 



B->1 



We find: 



(2/ 2 -i)! 



(2M-J) 1 



^2 M -2Z 2 
24(i 2 +i3H Ha) 



2(4^-47 2 ) 



= 1 



A=0 



E a (V2) 2A Qt x Ql 



4A/-)4(7 2 -A) 



(A.49) 



2 A 

E 

J'=0 



2A \ \q2I 2 -2X 

J J 1 e \ e 2A -J+ 1 /2 



£-►1 



1?— »1 



1 



6 



6 • • '£a-i 



^ M -2/ m+1 ^Q 2 Ql AQ3Q2 



2/t-i' 



(aM£ m ->l) 



The derivative over variable e and inverse derivative over variable d are independent on the summation index j and 
we can transfer them out of the sum. The summation over index j gives: 



£:(? 

j=0 v 



d7 x J 



eU 



2 A 



(A.50) 
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where 



U 



6 



1 



AQ 2 Q 1 AQ 3 Q 2 
For summation over index j in (|A.49[) we obtain: 



£2 • • • £a-i 
AQaQa-iJ \AQ0Q1 



^-2i m+1 \AQ2Q1 AQ 3 Q 2 



+ 



£2 ■ ■ ■ £a-i 
^4QaQa-i 



2m' 



(A.51) 



2 A 



aiJSm ->l,e->l,#— ►! 



We stress that the object J7 is proportional to A. In the above relation we proceed by evaluation of the derivative 
over variable e. We have: 



Q2I 2 - 



2 A 



1 



£ 2A+l/2 



eU 



2A" 



E--1 



i=0 

With help of identity: 



min (2A,2/ 2 -2A) . . ,„ , . , 



(A.52) 



(2A + 1/2) 



2J 2 -2A-i 



we converted (|A.52I) to the final form: 



(A.53) 



a2/ 2 -2A 



1 



-2A+1/2 



eJ7 



2A' 



e-»l 



2 A 



(A.54) 



MIN 



E (-1) 4 



i=0 



M7N\ /maa\ / 2/ 2 -1/2 
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[7 + i9 



where MI AT = min (2A, 2J 2 - 2A), and MAX = max (2A, 2J 2 - 2A). 

The sum over index i can be expressed by help of Gegenbauer orthogonal polynomials C" following the relation 
(see [H, vol. I, page 634): 



i=0 



n\ ( a - n+ 1/2 A /a 



( 2 ;)"(-I)"^ 2 "(VT^) 



(A.55) 



We insert (|A.54[) into the relation (|A.49[) . and using the identity: 



(1/2) 2A (27 2 - 2A)! ( 2 2 t 2 _ 1 2 / A 2 )=d/2) 2/2 



20 



we find: 



(V2) 2/2 n D i 



i 



2/i-2/ m+ i 



1 



AQ 2 Q 1 AQ 3 Q 2 



6 ■ • ■ £a-i 
^4QaQa-i 



(A.56) 



E 

A=0 



2J 2 - 1/2 
2MIN 



Qt x Qt {l2 ~ x) 



1 / 77 \ MIN 

d 02/*-a7 a ^2/ 2 -2A + M/ J V ) I (c/ + 1?) 2A-M7iV /^M c MAX-MIN+l/2 { ^Z#/U) 

We see that the term (1/2) determine the force of divergence. This means, that the leading term have this term 
as great as possible and this situation set if I 2 = /i. In this case all indices i m — 0. This means that the derivative 
operators have no effects in this case. The argument of the Gegenbauer polynomials are going to infinity as 

y/l/A. For A finite, but A — > we can use the relation (see Szego [|J) 



lim x~ n C*{x) = 2" 



Inserting this into above relation we have approximatively: 



n + A-1 



(1/2), 



1 



1 



1 



2/2 \AQ 2 Qi AQ3Q2 



AQaQ 



A-l 



2/i h 



EH 



A=0 



Qt x Qt il2 ~ x) 



^2^-2/2 ^2I 2 -2X+MIN^ I gj + ^2\-MIN 

Taking into account, that I 2 = \x, the inverse derivative over variable d have no effect and we red: 



(A.57) 



^2^-2/2 ^2I 2 -2X+MIN^ (jj + ^2X-MIN _^ ^2n-2X+MIN ([/ + ^)2A W / \ 

Taking into account that J7 ~ A we have: 



tf->i 



1 



aq a q 



A-l 



2// 



(Qo + Qt) 



For finite but A — > we approximatively find: 

{C(A) 2 ^} 



„ 2^(1/2) f :n " lirriX 
2^,2 M v / ; 2 p I ^7 



Finally, for the leading divergent terms for /i — > 00 of the asymptotic series for 5a in quasi-continuum relation we 
have: 



2/i 



(A.58) 
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